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Despite recent developments in management theory, maintaining a manufacturing schedule re¬ 
mains difficult because of production delays and fluctuations in demand and supply of materials. 
The response of manufacturing systems to such disruptions to dynamic behavior has been rarely 
studied. To capture these responses, we investigate a process that models the assembly of parts 
into end products. The complete assembly process is represented by a directed tree, where the 
smallest parts are injected at leaves and the end products are removed at the root. A discrete 
assembly process, represented by a node on the network, integrates parts, which are then sent to 
the next downstream node as a single part. The model exhibits some intriguing phenomena, includ¬ 
ing overstock cascade, phase transition in terms of demand and supply fluctuations, nonmonotonic 
distribution of stockout in the network, and the formation of a stockout path and stockout chains. 
Surprisingly, these rich phenomena result from only the nature of distributed assembly processes. 
From a physical perspective, these phenomena provide insight into delay dynamics and inventory 
distributions in large-scale manufacturing systems. 


I. INTRODUCTION 

The twentieth century saw industrialized societies de¬ 
velop with the support of organized modern manufac¬ 
turing, accompanied by a rapid increase in demand and 
consumption of goods. Scientific management of produc¬ 
tion systems dates back to the early twentieth century 
[l] . After Taylor’s pioneering work, innumerable stud¬ 
ies have been undertaken to control, optimize, and pre¬ 
dict production flow in factories Q. These studies 
have contributed to the design and management of man¬ 
ufacturing systems. However, the description of produc¬ 
tion flows by these theories, including the queuing theory, 
remains unsatisfactory.especially for complex, dynamic 
production systems [J, [5| . A pivotal factor impeding our 
understanding of production flows is the dynamic prop¬ 
erties of a complex production system as a “many-body 
system.” In general, constraints pertaining to the vol¬ 
ume of components and the finite capacity (buffer) for 
each job make the system dynamics complex. 

This paper proposes a physical approach—suitable for 
dealing with dynamic phenomena—to the complex sys¬ 
tem dynamics. We use a simple model that captures 
the essence of the dynamics and makes it visible. In the 
context of nonequilibrium statistical physics, the asym¬ 
metric simple exclusion process (ASEP) [fl-Q has been 
vigorously studied as the most archetypal model of par¬ 
ticle flow with the exclusion (blocking) effect, and non- 
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trivial behaviors specific to nonequilibrium systems have 
been discussed. In addition, the connection between the 
queuing theory and ASEP has been reported in recent 
studies 1,0. Thus, we find the intersection of studies 
on manufacturing systems and physics. Note that in the 
supply-chain management field, it is known that even a 
simple system can lead to highly complex behavior, in¬ 
cluding chaos [ll], [l^ and the so-called bullwhip effect 
|13h17I | , which has also been investigated from a physical 
perspective IlsHrD ]. 

Hopp and Spearman [22| attempted the idea of using 
a physical approach called factory physics to study pro¬ 
duction flow. They successfully systemized the funda¬ 
mental knowledge about production systems with sim¬ 
ple mathematics and static analyses. However, recent 
physical approaches using stability analyses fl8H2l| and 
the formulation presented in this paper use a dynamic 
treatment, and are thus beyond the scope of the factory 
physics concept. The main goal of the present study is to 
uncover and understand dynamic phenomena observed in 
complex manufacturing systems that are beyond intuitive 
prediction. For this purpose, simple models are used to 
identify the correlation between phenomena and causes. 
Different from recent similar attempts [l8l42T| . we use 
a full discrete model—space, time, and inventories are 
discretized—to realize this approach. This enables us to 
understand the system at the microscopic level, disre¬ 
garding relatively less important factors (e.g., input and 
output buffers at each job, and prediction and adaptation 
mechanisms), and to focus on the nature of “assembly.” 

In this paper, we focus on a set of jobs in a single man¬ 
ufacturing unit that constitutes a supply chain. Each job 
corresponding to materials’ processing does not predict 
its supply and demand, but follows the state of neighbor- 
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ing jobs. In the absence of prediction, the bullwhip effect 
does not occur. The primary goal of this study is to reveal 
how “assembly” processes on a large production line af¬ 
fect the overall system. The system is perturbed by three 
types of fluctuations: demand fluctuations, supply fluctu¬ 
ations, and removal of defective products. The assembly 
process involves merging two or more production flows, 
where material provisions from these streams synchro¬ 
nize (couple) [231 ]. This significantly increases complex¬ 
ity when the system is under perturbations, thus making 
rigorous analysis of such system very challenging. 

In previous studies, a set of assembly processes has 
been represented abstractly by a network of jobs, buffers 
(nodes), and topology of parts flow (links) 0, |23|426| . 
Here, we reformulate an assembly system as an interact¬ 
ing particle system Q and observe its dynamic aspects. 
The collective behavior of exclusive particles moving in a 
discrete network has received considerable attention from 
physicists [27l43T| . However, the effect of particle coales¬ 
cence (assembly) on these systems has not yet been fully 
understood. 

The reminder of the paper is organized as follows. The 
the focal model definition is given in the next section. 
To understand the model in detail, we first focus on the 
dynamics of a restricted parameter set (Sec. Mil) and then 
consider variations in the fixed parameters (Sec. IIVI) . 
Finally, we summarize the results and discuss the outlook 
of the study in Sec. |V] 


II. MODEL 

Consider a directed regular tree network whose inde¬ 
gree and outdegree are k and 1, respectively (Fig. [T]). 
Parts (particles) are transported along links, and at each 
node, k different parts are assembled, generating a part 
for the next node. A node has a buffer of size b for each 
incoming part; that is, a node can contain b parts of one 
type at the same time. Hence, the stock status at a node 
is described by a set of parts for the k buffers. The net¬ 
work has seN assembly stages, resulting in the need for 
k s raw materials. To investigate the statistical properties 
of large-scale production systems, we assume that k s is 
large. 

Basic rules of particle transportation in the system are 
as follows: at each time step, an assembly node sends 
a part to the next node if and only if all the required 
parts corresponding to the incoming links are stored in 
its buffers, and the buffer in the next node for the product 
the current node is creating is not full. Here, this node 
state is referred to as the “production state.” If even 
one part required in the production is absent, the node 
does not send a part. This node state is referred to as 
the “stockout state.” Even if the necessary materials 
are available, the node cannot send the product to the 
next node if the buffer of the next node is full. This 
state is referred to as the “demand-deficiency state.” In 
this paper, these three states comprise the “operational 


state.” Noted that, regardless of the state, each node can 
accept a part only if the corresponding buffer is not full. 
These procedures are performed in parallel at each node. 
Thus, the transition of a particular type of parts in a node 
in stage-cr (tr = 2, 3, • • • , s — 1) is summarized as follows. 
In the production state, the transition is —1 if the part 
is not provided by its upstream process; otherwise, it 
is ±0. In the stockout or demand-deficiency state, the 
transition is ±0 if the part is not provided by its upstream 
process; otherwise, it is +1. At stage-s (the network 
“leaves”), k parts are provided, each with probability 1 — 
£i n . The end product is removed at stage-1 (the “root”) 
with probability 1 — e out . The probabilities £i n £ [0,1) 
and e 0 ut £ [0,1) are the error rate of supply and demand, 
respectively. These rates are related to the intensity of 
their fluctuations. 

The model is later generalized by introducing a defect 
rate w ou t £ [0,1) reflecting the possibility of producing a 
defective product, which is immediately abandoned and 
not sent to the next node. When a node is in the pro¬ 
duction state, the produced part is sent to the next node 
with probability 1 — w ou t; otherwise, it is removed from 
the system as a defective part. 


III. FUNDAMENTAL PHENOMENA 

In this section, we restrict ourselves to considering the 
case of (w ou t = 0, b = 2, k = 2) as the most fundamental 
example. These restrictions will be relaxed later. In the 
absence of the possibility of producing defective parts, 
parts flow in each link is conserved. 

The case of b = 2 is important for the following rea¬ 
sons. First, this is the minimum size that allows a sta¬ 
tionary production flow. If each node can hold only one 
inventory item, products can be sent once per two steps 
as in ASEP. Second, having surplus stock in each pro¬ 
cess is costly. Generally, in the production management 
field, reducing stock has been encouraged |d2l fdi| . Given 
cost considerations and conventional domain recommen¬ 
dations, the buffer size we have set is ideal if supply fluc¬ 
tuations are negligible. Third, setting the buffer size to 
2 leads to a unique stationary state of production; that 
is, when demand and supply are steady, all nodes are 
in the production state with a single inventory item in 
each buffer. This state is robust for small fluctuations in 
supply, which will be discussed later. 


A. Overstock cascade 

First, we focus on the e ou t = 0 and small e; n cases 
to understand the fundamental system behavior. Fig¬ 
ure [2] illustrates the propagation of the stockout state 
and the generation of overstock (full buffers). For the 
stationary state, where all nodes are in the production 
state and their buffers have one part each, the delay of 
a single part is input. A stockout caused by supply flue- 


3 



FIG. 1. (a) Example of the assembly network (k = 2,s = 5). We define blocks-A and -B as the blocks surrounded by purple 
(left) and green (right) lines, respectively, (b) Local transition rules for b = 2 and k = 2. When the focal node (at right 
bottom) has both types of parts and the buffer in its downstream node is not full (bl and b2), it can send a part (production 
state). If even one required part is absent (b3), the node does not send a part (stockout state). When the node is not in the 
stockout state (that is, it is capable of producing a part) and the buffer in the next node is full (b4) it does not send a part 
(demand-deficiency state). 


tuation travels downstream (from a = s to 1). A node 
in the stockout state induced by a shortage of the fo¬ 
cal part stops its production, while the rest of the parts 
are provided to the node as scheduled, resulting in the 
buffers becoming full. For these excess inventories, up¬ 
stream processes stop production. This is observed as the 
emergence of the demand-deficiency state, which propa¬ 
gates upstream. This is in contrast with the stockout 
state because all inventories are increased by one at the 
demand-deficiency node in the next step. This leads to 
a demand-deficiency state at its upstream nodes. These 
local dynamics result in an “overstock cascade.” After 
these disruptions, each node recovers to the stationary 
production state in the presence of steady input and out¬ 
put. In addition, pair annihilation occurs when the stock¬ 
out state and the demand-deficiency state encounter each 
other. At this point, the demand-deficiency state spreads 
downward, with the exception of the link corresponding 
to the stockout. These states cancel each other out be¬ 
cause overstocking is avoided due to the shortage of in¬ 
coming parts. Here, the stockout state acts as a time 
buffer for overstocking. In addition, when two stockout 
states converge at the same node, they behave as a single 
state after the merger. Thus, these interactions decrease 
the number of stockout-state nodes. 


B. Production rate 

Here, we investigate the production rate q , defined as 
the expected output value of end products in one time 
step. Note that because each part is used once for one 
end product, the flow of parts in each link is identical to 
the production rate. 

Figure |3Ia) shows the production rate for various error 


rates. Here, we set s = 8. The most remarkable charac¬ 
teristic is the steep decay of q for e out = 0 and small £i n , 
which is asymptotic to the linear expression 

q- 1 - k s e in . (1) 

When the error rate ei n is sufficiently small, the effects of 
two different stockouts and their propagations do not in¬ 
teract, i.e., the occurrence probability of this interaction 
is negligible compared with e; n . Hence, a part shortage 
generated at leaves with rate k s e i n directly reduces the 
production rate. This expression is illustrated in the in¬ 
set in Fig. [3](a). Since k s is a large number, the produc¬ 
tion rate decreases significantly, even with a small error 
supply rate. After this steep drop, q slopes gently down¬ 
ward because the interaction between the stockout and 
demand-deficiency nodes is no longer negligible, pair an¬ 
nihilation and coalescence occur, and the supply errors 
do not directly decrease the flow. These characteristics 
stem from the effect of finite buffer capacity (i.e., par¬ 
ticle exclusion). When there is no limitation on buffer 
size, the production rate is expressed as q = 1 — £j n . (See 
Appendix ICl for its derivation.) 

When e out is nonzero, the production rate q has a 
plateau of q = 1 — e ou t) and the production flow is lim¬ 
ited by demand. In other words, in this region, nodes 
are almost fully occupied with inventories, and wait for 
the output that takes place with a probability of 1 — £ out . 
After this plateau, the flow becomes limited by supply 
errors near a certain e- ln . As such, the system has two 
different phases, the demand-limited phase (DP) and the 
supply-limited phase (SP), which are shown in Fig. [3Kb). 
In each phase, the production rate depends on only the 
corresponding error rate; and thus the curves in Fig. [3Ka) 
merge into the case of e out = 0. In DP, the bottleneck lim¬ 
iting the flow is clearly demand (output point of the end 
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FIG. 2. Overstock cascade. Blue, yellow, and red circles denote nodes in production, stockout, and demand-deficiency states, 
respectively. Blue arrows indicate parts flow at the time step; other arrows indicate absence of flow caused by stockout or 
demand-deficiency. We consider a situation where the supply and demand of parts at the boundaries are stationary (ei n = 
e out = 0). Here, we show the response of the system to a supply failure at the leftmost node in stage-4. The stockout 
propagates downstream, generating demand-deficiency for the different branches. In contrast, the demand-deficiency state 
propagates upstream. After 10 time steps, the uniform stationary state is recovered. 


product). However, the problem is not as simple in SP. 
There are many supply inputs and production paths, and 
the bottleneck that limits the flow contains a few stock¬ 
out nodes. Moreover, these limiting nodes move tempo¬ 
rally. Thus an analytical description of the flow is very 
challenging. 

C. Distribution of the stockout state 

For each stage, we study the probability of finding a 
node in the stockout state, p so ,a- Here, we do not fo¬ 
cus on the demand-deficiency state found with proba¬ 
bility pdd,a because a complementary relationship exists: 
Pdd.er = 1 — q — Pso,< 7 ■ Moreover, the production state 
can be found with probability q , irrespective of the stage 
number. These probabilities characterize the situations 
in each stage. Simulation results are shown in Fig. [4] 

For a small supply error e; n = 10 -5 , where interaction 
among delay propagations is negligible, the distribution 
is monotonic and follows 

p so ,a oc [k(l - eout) 2 ]" 7 - (2) 

This can be understood by the following approximate 
evaluation, which ignores the probabilities of 0(ei n ). In 
DP with small ej n (<^C e out ), each node is in either state-1 
or state-2; all buffers have one part (state-1) or are fully 
occupied (state-2) because the occurrence of a stockout 
that breaks uniformity is nearly negligible. Note that 
these states are defined only by buffers of a focal node and 
are independent of that node’s operational state, which 
takes the next node’s state into account. After producing 
an assembled part, the node assumes state-1; i.e., state-1 
appears with probability q = 1 — e ou t- Hence state-2 can 
be found with probability 1 — q = e out . We then consider 
the propagation of the stockout state generated at leaves. 
The probability of finding the incoming stockout state to 


a node in stage-cr is expressed as kp so . a +iq(\ — e ou t) by 
ignoring the correlation between the nodes. The factor 
kp SOt<7+ i represents the probability of having an incoming 
stockout state, and q is the probability of finding state-1 
at the focal node. If its next node is in state-2 (which 
state is observed with the probability e ou t), the propagat¬ 
ing stockout state vanishes. Thus, the factor q( 1 — e ou t) 
reflects the probability that the stockout state can sur¬ 
vive in the node. On the other hand, the probability of 
a stockout going out of the node is simply p SOjCr . Hence, 
by equating them, we obtain p SOjtT = fc(l - e ou t) 2 Pso,a+i- 
For example, for k = 2, e ou t = 1 — -j= — 0.2929 gives uni¬ 
form distribution of the stockout state. (See Fig. 0Ka) 
£out = 0.3.) In addition, when e ou t = 0, the nodes are 
nearly in state-1, and by ignoring 0(e? n ), the probability 
of a node entering the stockout state is given by fcp so ,< T +i, 
which is included in Eq. ©• These evaluations agree 
with the simulation results. Note that these arguments 
are not applicable to general values of b. 

The distribution shows an intriguing characteristic 
when the supply-error rate is large (ei n = 0.5). The stock¬ 
out probability p so has a local minimum around stage-5 
for cases in SP, while it is monotonic in DP. In SP, succes¬ 
sive nodes are observed in the stockout state (hereafter 
“stockout chain”). 

In the nodes providing parts to the stockout chain, 
the demand deficiency is induced at each time step. Be¬ 
cause the demand-deficiency state is generated upstream, 
a stockout state that appears in the network leaves van¬ 
ishes as it proceeds downstream. This is why p s0jfT in¬ 
creases from a = 6 to 8. In contrast, once the stock¬ 
out state reaches downstream stages, it becomes diffi¬ 
cult to eliminate and tends to stay as a stockout chain, 
thereby increasing its probability of remaining because 
the number of demand deficiencies generated in down¬ 
stream nodes decreases. In this phase, the demand defi¬ 
ciency generated by the output error rate e ou t at stage-1 
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FIG. 3. (a) Production rate for various error rates. The expression for plateaus observed in DP, q — 1 — tout, is illustrated 
with bold lines. The blue line in the inset corresponds to Eq. JTJ. (b) Phase diagram. When (e ou t) is relatively small, 
the production flow is restricted by the demand (supply), leading to DP (SP). The phase boundary is schematically indicated 
with a white dashed line, (c) Stockout-path probability, which is defined as the probability of the system having at least one 
stockout path. 


is less dominant than that generated by the supply er¬ 
ror e; n . For this reason, the stockout probability p so ,ir 
increases as the stage approaches completion, i.e., p SOjtT 
drops as cr increases. The relationship between e ou t and 
the amount of increase near the root (stage-1) can be 
understood by considering the following equation 

Q = (1 Pso,l)(l Cout)- (3) 

Since in SP, q is independent of e ou t; the stockout prob¬ 
ability must increase to p so ,i = 1 — g(l — e 0 ut) _1 at 
this point. Note that Eq. m is satisfied for general 
e, n , Cout) k, b, and s. Even if w out is nonzero, the modified 
equation, q = (1 - p so ,i)(l - e ou t)(l - w out ) holds. What 
is important in this figure is that for a fixed e- ln , the states 
in upstream processes are only slightly affected by e ou t 
and downstream processes in SP. Hence, a high stock¬ 
out probability for downstream nodes is not essential to 
ensure a high production rate; however, it is effective in 
reducing the stockout at leaves, for example, by extend¬ 
ing their buffers. This finding provides a new perspective 
for the buffer allocation problem (dbUioj . 


D. Stockout path and inventory imbalance 

With a large supply error, the production flow is dom¬ 
inated by the stockout state path (stockout chain) run¬ 
ning from the first to the final stage (hereafter “stock¬ 
out path”). When a certain part is not provided to a 
leaf node for a number of time steps, its downstream 
processes run out of these parts, while buffers for other 
parts in the system become full. Once the stockout path 
is formed, it tends to remain at the same location for 
the following reasons. While waiting for the focal part, 
inventories are accumulated in the other nodes, and as 
a result, when the required part is provided, the focal 
path consumes it immediately, once again exhausting the 
part. The probability of finding this stockout path in 




Stage number (b) Stage number 


FIG. 4. Stockout probability in each stage for (a) small 
supply-error rate tin = 1CF 5 (exponential behavior Eq. ([2]) ) 
and (b) large supply-error rate ei n = 0.5. (nonmonotonic re¬ 
lationship). 


the system (stockout path probability; SPP) is shown in 
Fig. He), where the stockout path becomes salient when 
ei n > 0.5 for relatively small e ou t, and ei n > 0.6 for large 
e ou t- Judging from the similarity between e ou t = 0.0 and 
0.5 in this region of ej n (SP), it is supposed that SPP 
is nearly independent of demand fluctuations. However, 
while the system is in DP and for sufficiently large e ou t, 
the increase in SPP is impeded by excessive inventories 
caused by small demand. 


We plotted the time series of buffer occupancy of the 
left and right halves of the system (Fig. |T]) in Fig. [5] 
For small supply fluctuations, two time-series transitions 
behave similarly; however, as the fluctuation increases, 
they desynchronize. Although the production rate is un¬ 
dermined through this decoupling, note that the average 
stock level remains largely unchanged. 
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FIG. 5. Time series of buffer occupancy in block-A (purple) and block-B (green) for e ou t = 0 and different values of £i n 
(0.2, 0.5, 0.8). As ei n increases, the two lines decouple. 


IV. EFFECTS OF NETWORK SIZE, BUFFER 
SIZE, AND PRODUCT INSPECTION 

Now, we focus on the dependence of variations in k. 6, s, 
and Wout on the system dynamics. 


A. Degree k and buffer size b 

Keeping the other parameters fixed as s = 8 and 
w 0 ut = 0, we discuss here the variations in degree k and 
buffer size b. Figure [HI a) shows the flux diagrams for 
k = 1,2, and 4 and 6 = 1,2, and 100 cases (and thus the 
center panel in each block is the same as in Fig. [3]). Note 
that each k = 1 case (shown for reference) corresponds 
to a chain of production processes and does not have as¬ 
sembly processes. When 6=1, the system coincides with 
ASEP. 

Generally, a too small buffer size (6 = 1) for mitigating 
exclusion (blocking) effects leads to the exacerbation of 
particle flux. A large difference in fluxes can be observed 
between 6 = 1 and 2 because 6 = 2 is the minimum buffer 
that allows stationary flow. The mitigation of the exclu¬ 
sion effect explains the sharpness of the phase transitions 
between SP and DP for larger 6. A large buffer is also 
effective for mitigating the effect of disruption in particle 
flow, and thus contributes to increasing the flow. 

Degree k is directly connected to the steepness of flux 
at £i„ = 0 [Fig. da)] as discussed in Sec. IIIIBI When 
the buffer size is large enough to absorb disruptions, the 
flux diagrams become independent of k: 

q = min {1 - e in , 1 - e out } (4) 

See Appendix [Cl for its derivation. 

Figure [7] displays the stockout probability p so ,a , and 
the particle density is defined as 


where r* denotes the number of particles in a buffer in 
a node in stage-tr. In the presence of system symmetry, 
this value is independent of the buffer position. As s 


increases, the particle density increases, absorbing fluc¬ 
tuation in inventory arrivals and thereby reducing the 
stockout probability. However, compared with the large 
number of inventories a node has to retain, the effect is 
not significant. This implies that a policy of increasing 
buffers is relatively ineffective for increasing the produc¬ 
tion rate. Rather, a focus on synchronization of the in¬ 
ventory level between nodes is more relevant to achieving 
cost-effective operations. 


B. Stage number s 

By investigating s = 1 cases, where assembly is con¬ 
ducted only once in the production process, and compar¬ 
ing the results with s = 8 cases, we discuss the effect of 
the stage number s. Figure [^b) shows flux diagrams for 
s = 1. Note that, because in these cases, the system is 
small enough, the flux can be calculated analytically (see 
Appendices [A] and Q3]). Stage number reduction directly 
leads to the total number of provision errors at supply 
boundaries, increasing the particle flux when the supply 
side is limiting the flow. On the other hand, when e out is 
more dominant, flux is deteriorated compared with s > 1 
cases, because buffers in each stage contribute to absorb¬ 
ing fluctuations that originate in supply errors, and thus 
plateaus of q = e ou t disappear when the buffer size and 
e 0 ut are small. 


C. Defect rate u> ou t 

In the presence of positive w ou t, particle flow is jeop¬ 
ardized s times by defects generated in the production 
process. Hence, when particles flow smoothly, the flux 
is reduced to (1 — w ou t) s times. Figure [8] shows relative 
differences from this ratio, defined as 

5=£ ^r [(1 - w out rr\ (6) 

where ( 7 ( 0 ) denotes the production rate when the defect 
rate is zero. These figures indicate that increasing k (6) 
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FIG. 6. Flux diagrams for various degrees k and buffer sizes 6. The defect rate is set to Wout = 0. For b = 100, theoretical 
predictions are given by Eq. 0, while for 6=1 and 6 = 2, analytical expressions are derived in Appendices [A] and [Bl 
respectively. 
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FIG. 7. (a) Stockout probability and (b) density distribution for various degrees k and buffer sizes 6. The other parameters 
are fixed as ei n = 0.5, s = 8, and uw = 0. 


strengthens (weakens) the effect of defects. This is be¬ 
cause large k increases the randomness of parts provision 
to downstream nodes, which decreases effective produc¬ 
tion rates, while large b contributes to mitigating this 
randomness, especially when demand is restricted (i.e., 
a sufficient number of parts are stocked in downstream 
nodes). Interestingly, these effects are highly dependent 
on w ou t- When w 0 ut is small, particles are not often re¬ 


moved, and the system is demand-driven. In this state, 
inventories are well stocked in downstream nodes, rein¬ 
forcing the effect of buffers. In contrast, when iL> ou t is 
large, many particles are removed before reaching the 
end of the tree, and the buffer effect is no longer present. 
This explains the non-monotonic behavior of 5 in Fig. [5] 

Figure 0 illustrates the stockout probability and the 
density distribution for w ou t = 0.15. This is a sufficiently 














































































large value for the particle removal process to determine 
the system dynamics. Hence, when k > 2, these distri¬ 
butions are almost independent of e out and determined 
mainly by network structure and buffer size. Note that 
for a large buffer size, they coincide with the e ou t = 0 
cases in Fig. 0 which is consistent with the fact that 
fluctuations caused by particle removal can be absorbed 
by large buffers. The effect of particle removal is more 
significant in downstream stages (small a) since particles 
become scarce as they approach completion. For smaller 
values of ui out , the aforementioned characteristics do not 
change drastically between Figs. [7] and [5] 


V. DISCUSSION 

We modelled the transportation of parts in a factory, 
and revealed fundamental characteristics of production 
flow that are caused by the effects of exclusion and coa¬ 
lescence of particles under supply, demand, and particle 
removal fluctuations. 

First, we confirmed the propagation of a stockout and 
demand deficiency caused by a single supply error (Fig. 
H). With this knowledge, we then analyzed the charac¬ 
teristics of the end product production rate and found 
steep decay for small supply-error rates (Fig. [11(a)). In 
general, the system behavior can be categorized into two 
distinctive phases depending on its output and input er¬ 
ror rates: DP and SP (Fig. |U(b)). 

In SP, production flow is restricted by a moving bottle¬ 
neck, i.e., successive stockout nodes (the stockout chain). 
For larger error rates, the production rate is mostly con¬ 
trolled by the stockout path, i.e., the stockout chain pen¬ 
etrates all stages, which tends to stay at the same loca¬ 
tion. Corresponding to this stockout path, the inventory 
level in blocks of the production network largely fluctu¬ 
ates (Fig. [5J. This inventory imbalance can be eased by 
the control policy. For example, stopping production if a 
node is producing more parts than the current bottleneck 
process could mitigate symmetry breaking. This would 
reduce cost for holding inventories and can stabilize the 
system. Also note that this prediction is consistent with 
the spirit of the kanban production system (4ll - [43j . in 
which upstream processes are permitted to produce only 
when they receive orders from downstream processes. 

In addition, for input errors that are not small -the 
order of “small” e; n is characterized as the absence of 
interactions between two input errors— nontrivial dis¬ 
tributions of stockout states can be observed (Fig. (4[), 
which have a minimum value of the stockout probability 
in the bulk. At the boundaries, we find local maximum 
values, and one of these values at stage-s is essential for 
controlling the production rate. From these observations, 
it can be concluded that in SP, for production systems 
with a large supply-error rate, it is more cost effective to 
increase the buffer size at the leaves than at downstream 
nodes in the assembly network. 

Through simulations, we found that the effects of 


buffer size, 6; degree of tree, k ; and stage number, s, 
can be summarized as follows (Figs. [G] and [Sj: (i) In¬ 
creasing k reduces particle flux, which can be absorbed 
by a sufficiently large b. However, compared with costs 
for maintaining large buffers, the effect is not significant; 
that is, the effect of flux fluctuations is not effectively ad¬ 
dressed by solely increasing buffers (as long as it is larger 
than the minimum size that allows steady flow (b > 2)). 

(ii) a small stage number s leads to the deterioration of 
flux in DP, while it contributes to increasing flux in SP. 

(iii) In the presence of positive defect rate w out , large 
b is effective in increasing DP flux, while k should be 
smaller when the system is not congested by particles. 
These findings are worth considering as fundamental to 
improving understanding and designing the dynamics of 
real manufacturing systems. 

When fluctuations are present, the stock level in each 
production block does not synchronize, which may rein¬ 
force production system instability. From a macroscopic 
view, among a group of factories (companies) that pro¬ 
vide parts for end products, the bullwhip effect can be 
intensified by this uneven distribution of stocks, and vice 
versa. Note that the network-induced nonuniformity and 
bullwhip effects are substantially different phenomena; 
the former is characterized as horizontal inventory im¬ 
balance and the latter is characterized as a vertical in¬ 
ventory imbalance that appears only when the control of 
orders is based on the prediction of future demand. 

This paper is devoted to finding the fundamental prop¬ 
erties of the process by setting a regular homogeneous 
network. However, it is also interesting to investigate 
the effects of network structure on production flow, het¬ 
erogeneity of processing time, and buffer size, which will 
be the focus of future work. Furthermore, it would be 
interesting to consider the technology diffusion in sup¬ 
ply chain networks in the model J44(. Although these 
changes influence the system dynamics, the presented in¬ 
terplay among various types of fluctuations, number of 
neighboring nodes, buffers are expected to be observed 
in general. 

Although some problems remain, we believe that fur¬ 
ther study of the physical approach presented in this pa¬ 
per is promising and will facilitate effective understand¬ 
ing and management of production systems (e.g., under¬ 
standing of how to place buffers effectively and character¬ 
istics of the response to environmental variations could 
be useful for addressing actual problems including the 
scheduling problem and cost optimization in factories). 
In this concept, the next goal is to show comparison with 
empirical data and application of statistical physics to 
the system. 
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Appendix A: Production rate for (s = 1, b = 1). 


When end products are manufactured by a single assembly (s = 1) of k parts and each buffer can contain one 

part at most at the same time (b = 1), the production rate can be calculated thorough the average number of time 

steps required to accept k different parts, T. Since parts provision to each buffer follows an identical and independent 
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Poisson distribution, one obtains 


p(T = t ) = P{T < t ) - P{T < t - 1) 

= (1 — e T ) fc — (1 — e T_1 ) fc , (Al) 

OO 

T = ^r[( l-c-) fc -(l-c T - 1 ) fc ]. (A2) 

T—0 


Considering time for production, parts for an end product are output every [T + (1 — e 0 ut) 1 ] time steps on average. 
Thus, the production rate is expressed as 


For example, in Fig. [G] (b), we use 


qk = 


T 


1 ^out 


(1 Ci n )(l e ou t) 

2 Ci n tout 

(1 — e in)(l + Cin)(l ~ Gout) 

2 Cout “b 2ej n (l Cout) C| n 
1 

l+5e ln +3e? n + 10e? n +2ef n +4ef n 1 

l + On + ^-^-dn-^n ^ l^ont 


(A3) 


(A4) 


Appendix B: Probability distributions for s = 1, k = (1,2), and b = (1,2) cases. 

When the number of states realized in a system is sufficiently small, we can calculate their probability distributions 
exactly. Let P Tl ... Tk be the probability of finding a state {ri, • • • , rjc}, where r, indicates the number of particles in the 
ith buffer. With a probability vector P = (P 0 ... 0 , • • • ,Pb---b) T and a transition matrix T( b+ i)k x ( b +i) k , the probability 
distribution is obtained through solving the following equation: 

T(b+i) k x(b+i) k P = P- (Bl) 

Using the obtained distribution, we can calculate particle flux as 

9 = (1 — Cout) Pri---T k i (B2) 

rii<i<fe Tijto 

where the sum is taken over all possible configurations of particles, where at least one particle is stored in all the 
buffers. 

In the following equations, explicit expressions of this equation for each (ft, b) are shown. 


1 . fc = 1,6 = 1 


Eqs. m 

feS) (*)-(*) < B3 > 

and (IB2I) yield 

q = (1 - e out )Pr = (1 ~ : *X ] - c °"t) , (B4) 

" £in ^out 

which coincides with Eq. (IATI) . Note that we use abbreviations ejif = 1 — e; n and = 1 — e out in the transition 
matrix. 
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2. k = l,b = 2 


Eqs. (1 bTT) 


and (1B2I) yield 


Em CinCout 0 
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(B5) 


(B6) 


Eqs. (IBTI) 


3. A: = 2,6 = 1 
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and (1B2I) yield 


(1 + £in)(l — e in)(l — E ou t) 

) O,- r \ ^2 
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(B8) 


4 . k — 2,b = 2 


Eq. 


is written as 
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Though particle flux can be obtained from 

q = P\i + P 12 + P 21 + P 22 , 

its explicit expression is too complicated to be written here. 


(B9) 


(BIO) 


Appendix C: Production rate for b 


We can also calculate the production rate when the buffer capacity is infinite. Differing from the previous section, we 
here focus on an arbitrarily long time T' and count the number of parts provided to the boundary in that period, since 
the parts’ inputs are never blocked. When T' is sufficiently large, the expected number of end products manufactured, 
qT ', coincides with the number of the part that is provided fewest times. Thus, we evaluate the minimum number of 
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parts provisions to bottom buffers. We define Si to be the number of parts provided to buffer i (i = 1, ■ • • , k s ) in a 
time interval T', and S m j n = mini< i </. s {Si}. 


P(S min > m) = P(Si > m) 


2=1 


- 


m —0 


T 


m 


‘T—m 


(1 - emf 


(Cl) 

(C2) 


q = Y’ mP (‘ S 'min = rn) 

m =0 

l r 

= 7p'Yh m >m)~ P(S min >m + 1)] 

m =0 

m=l 
— 1 Qn 7 


(C3) 

(C4) 

(C5) 

(C6) 

(C7) 

(C8) 


where U(x) is a unit step function that returns 0 if x < 0 and 1 otherwise. Approximations (IC6I) - (IC8I) are exact when 
T' —> cso. The obtained result q = 1 — ei n is consistent with the fact that the production flow is totally controlled by 
a single stockout path, to which a new material is provided with a rate 1 — 
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